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The scaling behaviors of the human electroencephalogram 
(EEG) time series are studied using detrended fluctuation 
analysis. Two scaling regions are found in nearly every chan- 
nel for all subjects examined. The scatter plot of the scaling 
exponents for all channels (up to 129) reveals the complicated 
structure of a subject's brain activity. Moment analyses are 
performed to extract the gross features of all the scaling ex- 
ponents, and another universal scaling behavior is identified. 
A one-parameter description is found to characterize the fluc- 
tuation properties of the nonlinear behaviors of the brain dy- 
namics. 

PACS numbers: 05.45.Tp, 87.19.La, 87.90. +y 



In the study of electrical activities of the brain recorded 
by electroencephalogram (EEG) various methods have 
been used to extract different aspects of the neuronal 
dynamics from the scalp potentials. They range from 
the traditional linear analysis that involves frequency de- 
composition, topographic mapping, etc. 0J^], to time- 
frequency analysis that uses wavelet transform ||, to 
nonlinear analysis that is particularly suitable for learn- 
ing about the chaotic behavior of the brain or 
for quantifying physiological conditions in nonlinear- 
dynamics terms |^|,^| . In this paper we discuss the scaling 
behavior of the fluctuations in EEG in nonlinear analysis 
and show the existence of new features in brain dynam- 
ics hitherto unrecognized. Moreover, we propose a global 
measure of the spatio-temporal signals that has potential 
utility in clinical and cognitive diagnostics. 

Apart from being more suitable to analyze non- 
stationary time series, the study of scaling behavior em- 
phasizes the relationship across time scales and provides 
a different description of the time series than the con- 
ventional Fourier power spectrum. Conveniently it also 
liberates our result from the dependence on the magni- 
tude of the voltage recorded by each probe. We aim to 
find what is universal among all channels as well as what 
varies among them. The former is obviously important 
by virtue of its universality for a given subject; how that 
universal quantity varies from subject to subject would 
clearly be interesting. The latter, which is a measure 
that varies from channel to channel, is perhaps even more 
interesting, since that quantity has brain anatomical cor- 
relates once the scalp potentials have been deconvolved 
to the cortical surface. 



Our procedure is to focus initially on one channel at a 
time. Thus it is a study of the local temporal behavior 
and the determination of a few parameters (scaling expo- 
nents) that effectively summarize the fluctuation proper- 
ties of the time series. The second phase of our procedure 
is to describe the global behavior of all channels and to 
arrive at one number that summarizes the variability of 
these temporal measures across the entire scalp surface. 
This dramatic data reduction necessarily trades detail for 
succinctness, but such reduction is exactly what is needed 
to allow easy discrimination between brain states. 

The specific method we use in the first phase is de- 
trended fluctuation analysis (DFA). This analysis is not 
new. It was proposed for the investigation of correlation 
properties in non-stationary time series and applied to 
the studies of heartbeat [Q and DNA nucleotides Q . It 
has also been applied to EEG ||, but with somewhat 
different emphases than those presented here. Since the 
analysis considers only the fluctuations from the local 
linear trends, it is insensitive to spurious correlations in- 
troduced by non-stationary external trends. By examin- 
ing the scaling behavior one can learn about the nature 
of short-range and long-range correlations, which are a 
salient aspect of the brain dynamics from the viewpoint 
of complex systems theory. 

Let an EEG time series be denoted by y(t), where t 
is discrete time ranging from 1 to T . Divide the entire 
range of t to be investigated into B equal bins, discard- 
ing any remainder, so that each bin has k — noor(T/£>) 
time points. Within each bin, labeled b (b = 1, • • ■ , B), 
perform a least-square fit of y(t) by a straight line, y b (t), 
i.e., y b (t) = Linear-fit[y(t)] for (b— l)k < t < bk. That is 
the semi-local trend for the 6th bin. Combine y b (t) for 
all B bins and denote the B straight segments by 

B 

y(k, t) = Y,y b (t) 9(t - (b - l)k) 9(bk - t) (1) 

6=1 

for 1 < t < kB. Define 
1 



kB 



F 2 (k) = —J2ly(t)-y(k,t)f 



(2) 



F(k) is then the RMS fluctuation from the semi-local 
trends in B bins each having k time points, and is also 
a measure of the fluctuation in each bin averaged over B 



1 



bins. The study of the dependence of F(k) on the bin 
size k is the essence of DFA (tLB. If it is a power-law 
behavior 

F(k) oc k a , (3) 

then the scaling exponent a is an indicator of the power- 
law correlations of the fluctuations in EEG, and is inde- 
pendent of the magnitude of y(t) or any spurious trend 
externally introduced. 

Resting EEG data were collected for six subjects us- 
ing a 128-channel commercial EEG system, with scalp- 
electrode impedences ranging from 10 to 40 k£l. The 
acquisition rate is 250 points/sec with hardware filtering 
set between 0.1 and 100 Hz. After acquisition, T w 10s 
lengths of simultaneous time series in all channels are 
chosen, free of artifacts such as eye blink and head move- 
ments. At each time point, the data are re-referenced 
to the average over all electrodes. This approximates 
the potentials relative to infinity, and provides a more 
interpretable measure of local brain activity We in- 
vestigate the range of k from 3 to 500 in approximately 
equal steps of lnfc. 



and subjects, it is a feature of EEG that conveys an im- 
portant property of the brain activity and should not be 
ignored. 
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FIG. 2. F(k) vs k for the three channels in Fig. 1. The 
vertical scales of Ch 1 and Ch 2 are shifted upwards by 1.0 
and 0.5 units, respectively. 
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FIG. 1. A sample of EEG time series in three channels. 
The vertical scales of Ch 1 and Ch 2 are shifted upward by 
60 and 30 /iV, respectively. 

In Fig. 1 we show three typical time series y(t) in three 
widely separated channels for subject A, labeled 1-3, for 
brevity. While it is clear that both channels 2 and 3 
have substantial 10 Hz oscillations after 0.2 s, it is much 
less apparent that there exist any scaling behaviors in 
all three channels. The corresponding values of F(k) are 
shown in the log-log plot in Fig. 2. Evidently, the strik- 
ing feature is that there are two scaling regions with a 
discernible bend when the two slopes in the two regions 
are distinctly different. With rare exceptions this feature 
is found in all channels for all subjects. Admittedly, the 
extents of the scaling regions are not wide, so the behav- 
ior does not meet the qualification for scaling in large 
critical systems or in fractal geometrical objects. How- 
ever, since the behavior is so universal across channels 



To quantify the scaling behavior, we perform a linear 
fit in Region I for 1 < In k < 2.5 and denote the slope by 
ati, and similarly in Region II for 3.5 < Ink < 5.75 with 
slope denoted by a 2 - Knowing the two straight lines in 
each channel allows us to determine the location of their 
intercept, In k, which we define as the position of the bend 
in lnfc. We find that, whereas a\ and a 2 can fluctuate 
widely from channel to channel, k is limited to a narrow 
range in most subjects. The average values of In k for the 
six subjects are 3.45, 3.66, 3.10, 2.92, 2.60, and 3.12. 

Since the only quantity in our analysis that has a scale 
is fc, the size of the bin on the time axis in which fluctu- 
ations from the semi-local trend are calculated, the spe- 
cific scale k must correspond roughly to a particular fre- 
quency /. If the data acquisition rate is denoted by r, 
then / = r/n. For r = 250 points/sec, and ln/c = 3.45, 
we get / = 7.94 Hz. That is near the midpoint separating 
the traditional a (8-13 Hz) and 9 (5-8 Hz) EEG frequency 
bands. Region I is then at higher frequencies, Region II 
lower. Thus in this nonlinear analysis we have found the 
existence of a specific frequency in every channel that 
separates different scaling behaviors. A change in scal- 
ing exponent in physical systems is often attributable to 
distinct dynamical processes underlying the generation 
of the time series. An interesting question is how this 
finding relates to EEG neurophysiology. 

We now exhibit the values of a\ and a 2 of all channels 
for subject A in a scatter plot in Fig. 3. They vary in 
the ranges: 0.19 < a x < 1.44 and 0.018 < a 2 < 0.489. 
Whereas a± is widely distributed, a 2 is sharply peaked at 
0.1 and has a long tail. The value of a = 0.5 corresponds 
to random walk with no correlation among the various 
time points. For a =/= 0.5 there are correlations: Region 
I corresponds to short-range correlation, Region II long- 
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range, with k giving a quantitative demarkation between 
the two. It is natural to conjecture that high deviations of 
the a values from the average could be caused by patho- 
logical conditions. Since the two scaling regions have 
corresponding frequency bands, one would likely want to 
focus on Region II to study brain states characterized by 
marked S (1-4 Hz) and 9 (5-8 Hz) activity. Sleep and 
cerebral ischemic stroke are two such examples. 
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FIG. 3. Scatter plot of a.i vs ct\ for subject A. The three 
channels exhibited in Figs. 1 and 2 are shown as big circles. 

So far our consideration has focused on the temporal 
properties of the time series in each channel. The scat- 
ter plot of a 2 vs ol\ provides a view of the complicated 
activities in all channels over the entire scalp. While the 
detailed spatial structure may be of interest for relating 
the values of on to brain anatomy, there is also a need 
for a general, overall description of all the pairs on. A 
global measure for each subject could be of great use to 
specialists and non-specialists alike. To that end, we now 
consider the relationship among the on values, or more 
precisely, the fluctuation of aj. 

Let x be either cvi or a 2 , and N be the total number of 
channels whose on values are under consideration. Since 
no value of on has been found to exceed 1.5 in the subjects 
we have examined, we consider the interval < x < 
1.5. Divide that interval into M equal cells, which for 
definiteness we take to be M — 150 here. Let the cells 
be labeled by m = 1, ■ • • , M, each having the size 5x — 
1.5/M. Denote the number of channels whose x values 
are in the mth cell by n m . Define 



(4) 



It is the fraction of channels whose x values are in the 
range (to — 1) 5x < x < mSx. By definition, we have 
J2m=i = 1- In Fig. 4 we show as an illustration the 
two graphs of P m for subject A, three of whose EEG 
times series are shown in Fig. 1. The two graphs corre- 
spond to a\ and 012, and are, in essence, the projections 
of the scatter plot in Fig. 3 onto the a\ and «2 axes. 
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FIG. 4. The distributions P m for a\ and 02 for subject A. 
The bin size in a for this plot is 0.02. 

Instead of studying the complicated structures of the 
distributions P m themselves, it is more convenient to ex- 
amine the moments of P m . Thus we define the normal- 
ized moments 



M 
m — 1 



M 



i q p m / (]T mp m y, 



(5) 



m— 1 



where q is a positive integer, although it can be a con- 
tinuous variable. Since x can be either ct\ or a 2 , we 



shall use Pm^ and Pm' to denote the two distributions, 



,(2) 



and Gq 1 ^ and G)f', respectively, for their moments. Since 
these moments are averages of (m/m) q , where to is the 
average-m, they are not very sensitive to to itself. They 
contain the essence of the fluctuation properties of 0*1,2 
in all channels. 
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FIG. 5. The q dependence of In for subject A. 



In principle, it is possible to examine also the moments 
for q < 0, which would reveal the properties of P m at 
low values of to. However, the accuracy of our data is 
not too reliable for low-fc analysis, since the 60 Hz noise 
due to ambient electric and magnetic fields has not been 
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cleanly filtered out. In this paper, therefore, we restrict 
our study to only the positive q values. For high q, the 
large m/fn parts of Fm' 2 ' dominate G^' 2 K 

In Fig. 5 the g-dependences of In G q ' are shown for 
the distributions exhibited in Fig. 4 for < q < 10. They 
appear to depend on q quadratically. The relationship 
between G q X ^ and G q 2 ^ is, however, extremely simple, as 
can be shown by plotting them against each other in a 
log-log graph. In Fig. 6 we show the straightline behavior 
for three of the six subjects, all of whose EEG time series 
have been analyzed in the same method described here. 
Thus we can associate a slope i) to every subject, and 
conclude that there exists a universal scaling behavior 

Gf*(Gf)\ (6) 

This remarkable behavior is valid for all subjects exam- 
ined, but the exponent r\ varies from subject to subject. 
Thus we have discovered a measure that characterizes all 
the aj values of a subject. 
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FIG. 6. A log-log plot of G q 2 ^ vs for three subjects. 
The solid lines are linear fits of the large q parts of the plots. 

We label the six subjects whom we have examined as 
A-F. Their r\ values are given below behind the subject 
labels. 

A (2.82), B (1.22), C (1.74), D (7.68), E (6.39), F (3.29) 

At this stage of our investigation we have not yet ar- 
rived at a point where we can give a definitive correlation 
between the values of 77 and some specific aspect of the 
brain activity. It is, nevertheless, of interest to note that 
the subjects A-C were regarded as healthy control sub- 
jects, while D-F have each had recent encounter with 
ischemic stroke. 

Before any tentative correlations may be inferred, we 
emphasize the need for caution since many factors were 
involved in data acquisition, e.g., varying degrees of 
awakeness of the subjects. Clinical and cognitive appli- 
cation is not the main object of this paper. However, the 



intriguing possibilities of the method presented here sug- 
gest extensive application of the analysis to many more 
subjects and for further investigation on what the scal- 
ing behaviors of the spatio-temporal EEG reveal. Such 
behaviors can undoubtedly be useful in guiding future 
theoretical development. 

To summarize, we have studied the scaling properties 
of the fluctuations in the time series, and found a number 
of features previously unrecognized in power spectrum 
analysis. The characterization of the temporal behavior 
of each channel by just two parameters, ol\ and a 2 , has 
led to the discovery of the r\ exponent. Although r\ is ob- 
tained by considering the global property of all channels, 
we envisage its extension to the study of local properties 
by the application of the moment analysis to partitioned 
regions of the scalp. Future possibilities of this approach 
to the investigation of EEG signals seem bountiful. 
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